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Abstract. We present our measurements of the deuterium to hydrogen ratio (D/H) in 
QSO absorption systems, which give D/H — 3.40 ±0.25 x 10~^ based on analysis of four 
independent systems. We discuss the properties of two systems which provide the strongest 
constraints on D/H. We outline the systematic effects involved in measurements of D/H 
and introduce a sophisticated method of analysis which properly accounts for these effects. 
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1. Introduction 

The status of standard big bang nucleosynthesis (BBN) has been extensively 
reviewed in these proceedings and in the recent literature (Copi et al. 1995; 
Sarkar 1996; Hata et al. 1997; Cardah & Fuller 1996; Fuller & Cardall 1996; 
Schramm &: Turner 1997). It has been emphasized that a determination 
of the primordial deuterium abundance has the potential to constrain the 
predictions of BBN, test the standard model, and give a precise measure 
of the present-day baryon density. Recent measurements of D/H towards 
bright, distant QSOs are now realizing this potential, and are already giving 
improved constraints on the models of BBN. (Tytler &: Buries 1997). 

Hogan (1998) gives a thorough overview of the present status of extra- 
galactic D/H measurements; here we focus on measurements of the systems 
for which we have obtained high-quality spectra with the HIRES spectro- 
graph (Vogt et al. 1994) and Keckl 10-m telescope. 



2. The Measurements 



We have analyzed four high-redshift QSO absorption systems which place 
useful constraints on D/H. We show the limits on D/H in each of these sys- 
tems in Fig. 1. Two of the systems, towards Q125H-3644 and Q1759+7539 
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Figure 1. Constraints on D/H from four independent measurements in four separate QSO 
absorption systems. Shown are the most likely values and 68% confidence limits at the 
absorption redshift of the D/H systems. 

have larger uncertainties for different reasons. Q125H-3644 is a faint QSO 
(V=19), and over 10 hours of observing time yields a spectrum with modest 
signal-to-noise ratio (SNR = 10) at high-resolution. The SNR drops quickly 
at lower wavelengths, and we cannot extract useful information from the 
highest order Lyman lines as a result. The system towards Q1759+7539 is 
a complex system of absorbers with a very high neutral hydrogen column 
density, log N(H I) > 10-*^^ (all column densities expressed in cm~^ ). The 
analysis of that system gives only an upper limit on D/H, at the 95% con- 
fidence level; the constraints are obtained from comparing the line profiles 
of Ly/3 and Ly7 to the profile of Lya . The method used to extract the 
limits on D/H will be described below. These results are preliminary, and 
improved constraints will require more data to allow a robust analysis. 

Tytler et al. (1996) made the first measurement of low D/H in the 
absorption system at z = 3.572 towards Q1937-1009. We analyzed the high- 
resolution spectrum (8 hrs of exposure), which resolved the entire Lyman 
series up to Ly-19, as well as associated metal lines. By profile fitting the 
Lyman lines, with the position of the velocity components given by the met- 
al lines, we find D/H = 2.3 ±0.3 ± 0.3 x 10~^ (statistical and systematic 
errors). The largest uncertainty in the measurement is the neutral hydro- 
gen column density, log N(H I) = 17.94 it 0.06 it 0.05, and the uncertainty 
stems from the saturated Lyman profiles (discussed in detail below) . We then 
obtained a high quality low-resolution spectra from Keck with LRIS (Okc 
et al. 1995), which gave better sensitivity shortward of the Lyman limit, to 
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Figure 2. Lyman limit region of Q1009+2956 at z = 2.504. 

Top: HIRES spectrum (FWHM — 8 km/s) showing Lyman lines 11 through 24, and model 
fit with three components and total H I column density, log N(H I) — 17.31. 
Bottom: Lick spectrum (FWIfM = 4 A) showing residual flux and model fit shortward of 
Lyman limit. The solid line shows the best fit model with log N(H I) — 17.39, and dotted 
lines represent 68% confidence levels. Absorption shortward of 3135 A is due to a lower 
redshift Lyman limit system. 
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directly measure the total N(H I) in the system and therefore place better 
constraints on D/H. Utilizing both the high and low-resolution spectra, we 
find log N(H I) = 17.86 it 0.02 by a direct measurement of the optical depth 
shortward of the Lyman limit at 4200 A (Buries k Tytler 1997a). With this 
constraint and a more sophisticated fitting procedure, we measure D/H = 
3.3 ±0.3 X 10-5 (Buries & Tytler 1997b). 

We discovered an absorption system at z = 2.504 towards Q 1009+2956 
ideal for a measurement of D /H. This system has a lower hydrogen column 
density, log N(H I) = 17.39 ± 0.06. The highest order Lyman lines become 
unsaturated, which yields a precise measurement of N(H I) in both low 
and high resolution spectra (Fig. 2). Over twelve hours of Keck+HIRES 
produced a very high quality spectrum of the entire Lyman series, resolving 
the entire series up to Ly-22. We find strong evidence for contamination of 
the deuterium Lya; absorption feature, which unfortunately introduces the 
largest uncertainty in the measurement. With the contamination included 
in the analysis, we measure D/H = 4.0 ± 0.7 x 10"^ (Buries & Tytler 1997c). 

The four independent systems support a low primordial abundance of 
deuterium, and together give D/H = 3.4 it 0.25 x 10"^. If this represents the 
primordial value, nucleosynthesis calculations from standard BBN models 
with three light neutrinos give the baryon-to-photon ratio, 77 = 5.1 =b 0.3 x 
lO"^'' and the baryon density, r^fe/i^QQ = 0.019 ± 0.001. We arrive at very 
small statistical errors (10% at 95% confidence), so we must confront the 
existence of systematic effects. In the remainder of this paper we present 
a sophisticated method that we use in order to properly account for the 
systematic effects currently known. 



3. An Improved Method to Measure D/H 

Many groups have reported detection of deuterium in absorption along the 
line of sight to distant QSOs (Songaila et al. 1994; Carswell et al. 1994; 
Tytler et al. 1996, Wampler et al. 1996; Rugers & Hogan 1996a,b; Webb et al. 
1997; Buries & Tytler 1997b,c). Here, we put forward an improved method 
to measure D/H, which has many advantages over the current, widely used, 
methods. 

The goal is to extract from the spectrum of a QSO absorption system 
the most likely value of D/H, and the confidence intervals about that value. 
We are still required to model the D/H absorption system (DHAS) with 
a finite set of parameters (based on physical assumptions) to compare to 
the observed spectrum. With a given model, we can quickly minimize the 
difference between model and data, parameterized as x^, and use the final 
parameters to find the most likely value of D/H. But the uncertainty in D/H 
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Figure 3. x functions of D/H in two separate absorption systems. The Relative x is 
shown to facilitate a direct comparison. The minima represent the most likely values of 
D/H in each system, and Ax^ = 4.0 represents 95% confidence limits. 



depends on correlations in the parameters, and is non-trivial to extract from 
the best fit solution. 

The function depends on the N parameters included to model the 
observed spectrum. Instead of mapping the function in the full N dimen- 
sional parameter space, we choose to map x^ versus a single parameter, 
D/H (known as the "ridge" method). Towards Q1937-1009, for example, we 
compare a model with 234 free parameters to regions of the high-resolution 
spectrum which total 1298 pixels. We note the inventory of free parameters: 
3 main components in Lyman limit system with 4 parameters each (N(H I), 
^ahsi T, htur) + 9 Spectral regions with a total of 29 parameters for the con- 
tinuum -I- 64 other H I lines with 3 parameters each (N(H I), Zahs-, hotal) 
which are near and affect the absorption lines of the D /H system + 1 param- 
eter for D/H: 3x4 + 29 + 64 x3-1-1 = 234. To construct the minimized 
X^ function versus D/H alone, we choose a value for D/H, fix D/H for the 
iteration, and allow the other 233 parameters to freely change to minimize 
X^- We choose and fix a new value for D/H, and minimize x^ again. The 
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Figure 4- Top: Lya region of z — 2.504 system towards Q1009+2956. The histogram shows 
the pixels of the HIRES spectrum normalized to the initial continuum estimate (solid line 
at unity). The la errors per pixel is shown as the solid line near zero. The smooth dashed 
line near unity is the best fit continuum with five degrees of freedom. The model profile is 
composed of the main H I (dot-dashed) and D I (dashed), and contaminating H I (dotted). 
Bottom: Lya region of 2: = 3.572 system towards Q1937-1009. The dashed line shows the 
best fit continuum level. 
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results for both Q1937-1009 (this example) and Q1009+2956 are shown in 
Fig. 3. In Fig. 4 (bottom panel), note the large number (19) of H I lines 
which must be included in the Lya region of Q1937-1009. 

In this paper we present the difficulties inherent to measurements of D/H 
in QSO absorption systems, and discuss techniques we have implemented to 
overcome the systematic effects. We outline an improved method to measure 
D/H, which accounts for the total uncertainty, in contrast to the statistically 
uncorrelated errors of D and H column densities summed in quadrature. We 
show that the best two current measurements of D/H are complementary 
and give a robust measure of the primordial value, D/H = 3.4 it 0.25 x 10~^. 

4. Systematic uncertainties in D/H measurements 

In this section, we describe the systematic uncertainties which are known to 
exist in measurements of D/H in QSO absorption systems. We explain how 
we account for the systematics, with the assumption that the nature of the 
effect is understood. 

4.1. Intrinsic QSO continuum 

The observed QSO spectra are combinations of the intrinsic spectrum of the 
source with intervening absorption features. In the measurement of D/H, we 
are solely interested in the intervening features. We attempt and require the 
removal of the intrinsic QSO spectrum (continuum) to analyze the absorp- 
tion features in the observed spectrum. The QSO continuum changes on 
much longer wavelength scales than the narrow absorption lines. A polyno- 
mial of low-order traces the smooth variations of the QSO continuum and 
gives the normalization of the observed spectrum to the QSO continuum. 
This procedure (called "continuum fitting") is not exact, and even worse, not 
well defined. However, abundance measurements in QSO absorption systems 
have not taken into account the uncertainties introduced by subjectively fit- 
ting a fixed continuum to the observed spectrum. Once a continuum fit is 
complete, the solution is assumed, for all practical purposes, correct with 
no uncertainty. We have overcome this problem by assuming the continuum 
is not fixed or well determined. As a part of the abundance measurements, 
we include free parameters to account for the unknown shape and level of 
the continuum in regions of interest in the observed spectrum. We still fit 
an initial continuum (a first guess), but this is not assumed absolute or 
correct, only a likely estimate. The procedure has been outlined in Buries 
& Tytler (1997b), and should be adopted for all abundance measurements 
which could be influenced by continuum uncertainties. 

As examples, we show the Lya regions of Q1009-I-2956 (top) and Q1937- 
1009 (bottom) in Fig 4. The initial continuum estimate is the solid line 



tytler.tex; 5/02/2008; 17:27; no v.; p. 7 



8 



SCOTT BURLES AND DAVID TYTLER 



at unity and is used to normalize the spectrum before profile fitting. The 
smooth dashed line near unity is the resulting best fit continuum with five 
free parameters. The level of the final continuum lies near the initial estimate 
on the D/H Lya feature, which confirms that the initial estimate was fairly 
well chosen. Allowing for a free continuum does not necessarily change the 
central value, but ensures that the final uncertainties take properly into 
account the uncertain nature of the intrinsic QSO continuum. 

4.2. Hydrogen Contamination of Deuterium 

The absorption features of the Lyman series of H I and D I are modeled to 
measure D/H in QSO absorption systems. These features reside exclusively 
in the Lya forest and are subject to blending, and therefore, contamination 
from other Lyman absorption. The line density of the Lya forest rises steeply 
with redshift, so measurements at higher redshift, z > 2, are subject to more 
contamination from random, unrelated Lya lines. The a posteriori proba- 
bility of unrelated H I contamination can be estimated from the statistical 
distirbution of Lya forest lines. (Steigman 1994, Tytler et al. 1996, Tytler 
& Buries 1997, Jedamzik & Fuller 1997, Hogan 1998). On the other hand, 
structure in the absorption systems will increase with decreasing redshift, 
so contamination of measurements at lower redshift, z < 1, is dominated 
by Lyman absorption correlated on scales of the D-H velocity separation, 82 
km s""^ . We must test for the presence of contamination in all measurements 
and ensure that the final results properly account for its existence. Contam- 
ination always gives lower values of D/H. A measurement which does not 
include contamination in the analysis gives an upper limit on D/H. In Fig. 
4, we show the individual profiles of H I (dot-dashed) and D I (dashed) in 
the Lyman limit system and unrelated H I contamination redward of the 
D I profile (dotted line). Contamination was included in analysis towards 
both Q1937-1009 and Q1009+2956. 

4.3. Saturated Absorption features 

In measurements of D /H, one must select systems with high neutral hydro- 
gen column densities, log N(H I) > 17, to be sensitive to D/H values less 
than 10~^. In these systems, most of the lines of the Lyman series have 
saturated (optical depth, r > 3) line centers. We obtain spectra with high 
signal-to-noise ratio (SNR) of the region shortward of the Lyman limit to 
independently constrain N(H I). It is possible to obtain a direct measure- 
ment of the total H I column in systems with 17 < log N(H I) < 18, 0.6 
< tll < 6.0, where tll is the optical depth of Lyman continuum absorp- 
tion at 912 A (rest frame). The residual flux shortward of the Lyman limit 
provides a direct estimate of N(H I), independent of other characteristics of 
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the system (e.g. velocity dispersion or deuterium column density) . To obtain 
an accurate measurement, we must account for unrelated Lyman absorption 
from other systems and the uncertainty of the intrinsic QSO continuum near 
the Lyman limit (see sec. 4.1). The majority of the Lyman absorption fea- 
tures can be specified by their related absorption features in high-resolution 
spectra (at higher wavelengths). The remaining features are statistically 
included from well determined Lya forest distributions (Kirkman &; Tytler 
1997). The QSO continuum near the Lyman limit can be determined to 
better than 10% by utilizing high-resolution spectra (c.f. Buries & Tytler 
1997a). The uncertainty in the N(H I) measurement is proportional to the 
continuum uncertainty, A N(H I) = A (cont) /tli,. In Fig. 2, we show the 
Lyman limit region of the z = 2.504 D/H system in two separate spectra 
of Q1009+2956. The high resolution spectrum (top, FWHM = 8 km s"^ , 
0.09 A) shows the resolved high order Lyman lines up to Lyman-22. The 
Lyman lines become unsaturated at the highest orders, and give an optical 
depth in the continuum, tll = 1-3. The low resolution spectrum (bottom, 
FWHM = 4 A) has better sensitivity and shows significant flux shortward 
of the Lyman limit at 3200 A. 

A one-parameter model, N(H I), is compared with the low resolution spec- 
tra to find the most probable value and confidence intervals of N(H I) (Buries 
&: Tytler 1997a). The features seen in the model fit in Fig. 2 are higher- 
order Lyman lines (Ly/3 and above) associated with Lyman lines measured 
at longer wavelengths in the high resolution spectrum. Lya absorption with- 
in 1.57 < z < 1.64 is drawn from a distribution and included in the Monte 
Carlo analysis which gives the maximum likelihood of N(H I): log N(H I) 
= 17.39 ±0.06, which includes the uncertainties from extrapolation of the 
unabsorbed QSO continuum. 

4.4. The Case of Mesoturbulence 

Levshakov et al. (1998, and references therein) have studied the effects of 
correlations in turbulent velocity fields (mesoturbulence), on line formation 
in the Lya forest. The standard line profiles (Voigt) are, in fact, the spe- 
cial case of mesoturbulence when the correlation length goes to zero. A 
substantial amount of work has been published to understand the impact 
of mesoturbulence on the abundance measurements in absorption systems, 
with emphasis on D/H (Levshakov &: Kegel 1997, Levshakov et al. 1997, 
Levshakov et al. 1998). We would like to point out a few observational con- 
straints on mesoturbulence, and the limited role it plays (in our view) in 
measurements of D /H. 

(1) The turbulent component should be consistent with all lines in the 
absorption system. Mesoturbulent profiles will significantly differ from Voigt 
profiles only when the turbulent component is comparable or larger than 



tytler.tex; 5/02/2008; 17:27; no v.; p. 9 



10 



SCOTT BURLES AND DAVID TYTLER 



the thermal component. Both components contribute to the total veloci- 
ty dispersion (htotal) which is measured in high resolution spectra: bl^^^^ = 
^meas = ^turb + ^Iherm' Comparisons of the metal line widths (e.g. C, Si) to 
the hydrogen and deuterium line widths place strict limits on the turbulent 
component. The heavy metal lines (Si and above) are dominated by turbu- 
lence at temperatures T < 10^ K and their velocity dispersion usually gives 
a good measure of the turbulence in the system. In the two systems we have 
measured D/H, towards Q1937-1009 and Q1009+2956, the metal Hues are 
very narrow, which give very low velocity dispersions, 6(Si) < 6 km s^^ . 
The hydrogen lines are much wider, with typical widths fe(H) ~ 20 km s^^ 
, which shows that thermal, uncorrelated, line formation is dominant in the 
optically thick Lyman lines. 

(2) The line profiles of mesoturbulence can be distinguished from Voigt 
profiles by fitting the entire Lyman series (Levshakov &: Kegel 1996). The 
requirement for invoking mesoturbulence should be determined by the obser- 
vational data. If the Lyman series lines show evidence for mesoturbulent 
profiles, then one should take it into account. In D/H systems to date, sim- 
ple Voigt profiles can explain all of the profiles in the Lyman series and do 
not require the mesoturbulent model. 

(3) Mesoturbulence is likely to play an important role in a wide variety of 
absorption line studies. The hydrogen lines in QSO spectra which do show 
deuterium belong to a class of absorbers with the narrowest intrinsic widths, 
and are therefore subject to the least amount of turbulent broadening. On 
the other hand, optically thick absorption lines of heavy metals (for instance, 
in damped-Lya systems) are prime candidates for studies of mesoturbulence. 

4.5. Line blending 

The D/H absorption systems do not have simple, symmetric profiles which 
can be modeled as a single component. Analysis of the absorption profiles 
of the Lyman lines and associated metal lines requires multiple components. 
The components are usually separated by less than the intrinsic line widths 
of the individual components, which gives rise to blended profiles of the 
Lyman series. The severe blending represents a loss of information. The 
parameters describing individual components are poorly constrained, but 
those describing the entire system can be characterized and measured. Here 
lies the problem: the set of parameters describing the multiple components 
of the D/H system are not orthogonal when the components are blended. 
The uncertainties in the parameters are strongly correlated, and the variance 
in the parameters does not represent the total uncertainty. Therefore, one 
cannot measure the individual column densities of D and H accurately in 
each absorption component. 
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